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Abstract 

We study the glass transition of binary mixtures of dipolar particles in two dimensions within the framework of mode- 
coupling theory, focusing in particular on the influence of composition changes. In a first step, we demonstrate that 
the experimental system of Konig et al. [Eur. Phys. J. E 18, 287 (2005)] is well described by point dipoles through 
a comparison between the experimental partial structure factors and those from our Monte Carlo simulation. For 
such a mixture of point particles we show that there is always a plasticization effect, i.e. a stabilization of the liquid 
state due to mixing, in contrast to binary hard disks. We demonstrate that the predicted plasticization effect is in 
qualitative agreement with experimental results. Furthermore, also some general properties of the glass transition 
lines are discussed. 

Keywords: Glass transitions of specific systems, Theory and modeling of the glass transition, Colloids 
PACS: 64.70.P-, 64.70.Q-, 82.70.Dd 



1. Introduction 

The study of supercooled liquids close to the glass 
transition is a rapidly developing field in modern 
physics. For both experimental and theoretical investi- 
gations, colloidal dispersions are often chosen as model 
systems. Compared to atomic liquids, the main advan- 
tages of such model systems rely on their simplicity: the 
particle interactions can often be described theoretically 
by quite simple and well defined expressions, and ex- 
perimental investigations can be performed by applying 
standard optical methods in real space instead of neu- 
tron or X-ray scattering techniques. 

Supercooled liquids close to the glass transition ex- 
hibit complex relaxation patterns which do not follow a 
simple exponential decay in time. The most success- 
ful theory for the description of the dynamics of su- 
percooled liquids is the so-called mode-coupling theory 
(MCT) which has been developed and studied in great 
detail by Gotze and coworkers Jll for various model 
systems in three spatial dimensions (3D). The central 
prediction of MCT is a transition from a liquid into an 
ideal glassy state upon decreasing temperature or in- 
creasing particle density below or above some critical 
value, respectively. 

It is well-known that many physical phenomena like 
equilibrium phase transitions depend strongly on the 
spatial dimensionality d. Consequently, it appears to 



be essential to study the ^-dependence of the liquid- 
glass transition as well. As far we know, the first ex- 
perimental approach to this issue is the investigation of 
Konig et al. J2J] of a two-dimensional (2D) binary mix- 
ture of super-paramagnetic colloidal particles which are 
trapped on a water-air interface. They interact via re- 
pulsive dipole potentials. Here, the magnetic dipoles 
have been induced by an external magnetic field per- 
pendicular to the water interface which renders the re- 
pulsions isotropic in the interface plane. Their results 
for the self-intermediate scattering functions measured 
by video microscopy clearly exhibit the stretched re- 
laxation patterns of glass-forming liquids and the cor- 
responding strong slowing down of the dynamics when 
increasing the magnetic field (which plays the role of an 
inverse, effective temperature). This model system was 
also used for both experimental and theoretical studies 
on crystallization effects in the high-coupling limit, i.e. 
for large magnetic fields (see ^ and references therein). 
A further interesting phenomenon (found experimen- 
tally and corroborated theoretically) pertaining to the 
same model system in the liquid state is a partial clus- 
tering of the smaller particles in a sponge-like topology 

a. 

The first MCT calculations on a model system in 2D 
were performed by Bayer et al. ^ for monodisperse 
hard disks. An ideal glass transition has been found at 
a critical packing fraction ip c . This study was extended 
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to binary mixtures of hard disks [6] characterized by the 
total packing fraction ip, the concentration of the smaller 
disks x s , and the size ratio 6. The same four mixing ef- 
fects occur as have been reported before by Gotze and 
Voigtmann (01 for binary mixtures of hard spheres in 
3D: (i) for small size disparities the glass is stabilized, 

(ii) for large size disparities the liquid state is stabilized, 

(iii) upon increasing the concentration of the smaller 
particles the plateau values of the normalized correla- 
tion functions for intermediate times increase for almost 
all wave numbers accompanied by (iv) a slowing down 
of the initial part of the relaxation of the correlators of 
the bigger particles towards these plateaus. Further, it 
was shown that the glass transition diagram for binary 
hard disks, i.e. the critical packing fraction <p c as func- 
tion of x s for several 8, strongly resembles the corre- 
sponding random close packing diagram. This fact is an 
indication for the validity of the MCT approximations 
also in 2D. 

In the present study, we will apply MCT to a 2D bi- 
nary liquid of dipoles. Our aim is twofold. First, we 
want to explore how far the mixing effects found for 
hard spheres 10] and hard disks |@] have some generic 
features. Second, motivated by the experiments on 
quasi-2D binary mixtures of super-paramagnetic par- 
ticles J2], we want to compare our MCT results with 
available data. As a byproduct, we propose a simple 
empirical ansatz for the bridge functions j9'] appearing 
in the closure to the Ornstein-Zernike (OZ) equation 
which allows for a rather precise determination of the 
equilibrium structure. 

The organization of our contribution is as follows. 
Sect.|2]will contain (i) a justification of the model used 
to approach the experimental system, based on a com- 
parison of the experimental and theoretical results for 
the static structure factors, (ii) a concise review of MCT 
equations in d dimensions, and (iii) the discussion of 
some properties of the glass transition lines. Our results 
from MCT will be presented and discussed in Sect. [3] 
and the final Sect. [4] summarizes and gives some con- 
clusions. 



2. Model system and mode-coupling theory 

In this section, we will introduce the model being 
investigated, discuss some of its structural properties, 
and shortly review the corresponding MCT equations. 
For more details on the latter the reader may consult 
Ref. H. 

Let us summarize some basic notation. We make 
use of the compact mathematical notation introduced in 



Ref. Bold symbols A, B etc. denote arrays of ma- 
trices. Their components A*, B^ labeled by subscript 
Latin indices are M x M matrices. Their elements A k , 
B°f \ also denoted by (A)? 8 , (B)f when appropriate, are 
indicated by superscript Greek indices. Matrix prod- 
ucts are defined component-wise, e.g. C = AB means 
Ck - A/cBk for all k. A is called positive (semi-)definite, 
(A > 0) A > if this is true for all A k . 

2.1. The model system 

Our starting point is the experimental model system 
studied by Konig et al. Technical details for the 
experimental setup can be found in Ref. Jit]. The au- 
thors consider a binary mixture of super-paramagnetic 
colloidal particles trapped on a water-air interface of a 
pending droplet, so that particle motion is effectively 
only possible in 2D. By applying an external magnetic 
field perpendicular to the liquid interface, they can vary 
the strength of the repulsive dipolar interaction forces 
between the particles. Temperature T, total particle 
number density «, radii R a and susceptibilities Xa of 
the smaller (a = s) and bigger (a = b) particles are 
kept constant to good accuracy. To observe vitrifica- 
tion, the magnetic field is increased. Accordingly, the 
magnetic moments of the particles and the repulsions 
between them grow and in turn the correlations between 
the particles become stronger. 

The colloidal dispersion is strongly diluted. To be 
more precise, the average distance between the parti- 
cles exceeds the diameter of the bigger particles roughly 
by a factor of four H, which yields the estimate ip < 
(4 ■ 2R b )- 2 nR 2 h = n/64 * 0.05 for the total 2D-packing 
fraction on the interface plane. Therefore, collisions of 
the hard cores of two particles practically never occur, 
at least for low enough temperatures. Under this condi- 
tion, the thermodynamic equilibrium state of the experi- 
mental system is well determined by the induced dipole 
potentials, only. 

Idealizing this picture, we assume a binary mixture of 
point particles in 2D interacting via the induced dipole 
potentials 



u a \r) = 



4k r 5 



(1) 



jUo is the vacuum permeability and B denotes the applied 
external magnetic field. Although the hard core interac- 
tions have been neglected, the two species of particles 
shall be called "big" {a — b) and "small" (a = s) in the 
sense that^;, > x.i- 

The thermodynamic state of the considered system 
with fixed Xb and Xs depends on four variables: n, T, 
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the concentration of the smaller dipoles x s , and the mag- 
netic field B. For a monodisperse dipolar system with 
susceptibility xo it is well-known that the dependence 
on T, n, B is given by the dimensionless coupling pa- 
rameter 

iq =n — -r^r- (2) 



4n k B T 

Since we want to study mixing effects on the ideal glass 
transition, we have to compare e.g. the critical temper- 
ature T c (x s ) or equivalently the critical filed B c (x s ) for 
the binary system with concentration < x s < 1 with a 
corresponding monodisperse system. As we will argue 
now it does not make sense for all x s to simply choose 
for that comparison the monodisperse system with sus- 
ceptibility xb vtXs- 

The magnetization per site m is a relevant physical 
quantity. For the binary system it is mbi, wry {x s ) = [x^-l- 
(1 - x s )xb]B. Now we will consider a monodisperse 
system with a susceptibility x( x s) f° r oil point dipoles. 
Then it is m mono (x s ) = x( x s)B. For a comparison of 
the binary and the monodisperse system it is reasonable 
to require that both have the same magnetization. From 
this condition we getx(x s ) = x s Xs+(l-Xs)Xb an d finally 
the generalization of Eq. (0 to binary dipolar systems in 
2D: 



T = (7m) 



3/2 /"Q fcsXs + (1 ~ x s )xb] 2 B 2 



4/r 



k B T 



(3) 



which is the dimensionless interaction parameter chosen 
by Konig et al J2l- In the following we will use this pa- 
rameter, the susceptibility ratio 6 = XslXb ^ 1- an d the 
concentration x s of the smaller dipoles as independent 
control parameters. The prefactor 7T 3 / 2 in Eq. (0 was 
introduced to match the convention used in Ref. |2J]. 

2.2. Static structure functions 

The interaction potentials for the specified model sys- 
tem enter the MCT equations only via the static struc- 
ture factors. Thus, we need some accurate and efficient 
method for a systematic calculation of these objects. We 
have at least two possibilities. 

First, we can make use of the well established integral 
equation theory for the direct correlation functions c" k 
based on the OZ equation |@] for mixtures 



h — c + ncxh 

a/3 



(4) 



where x" k P - x a 8 a p and h°f are the total correlation func- 
tions for given wave number k. 6 aB denotes the Kro- 
necker delta. To obtain e.g. c, we need a closure rela- 
tion. The general ansatz fl9[] for this reads 
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Figure 1 : (Color online) Partial structure factors for a binary mixture 
of dipolar particles in 2D at T = 110, 5 as 0.1, and x s = 0.48 mea- 
sured experimentally by Ebert (T^l (solid lines) and the correspond- 
ing MC results (dashed lines) obtained with 600 particles without data 
smoothing. 



u eff {r) 



u a Kr) B 



(6) 



(r), 



(5) 



g a P(r) = ffP(r)+l are the pair distribution functions and 
d Q (r) are called the (still unknown) bridge functions. 
Every specific closure relation relies on some (uncon- 
trolled) approximation scheme for c?z(r). The advan- 
tage of using an integral equation theory is, however, 
that data with high numerical quality can be produced. 
There is no statistic noise and there are no finite size 
effects like in atomistic computer simulations. Further- 
more, crystallization effects can explicitly be excluded. 
Having found c, the static correlation matrix S follows 
from (S -1 )^ = 5 oj3 /x a - ncf . 

Second, we can perform an atomistic computer sim- 
ulation. Since we are only interested in static quantities, 
the most efficient choice is the Monte Carlo (MC) tech- 
nique. It yields exact numerical results in the sense that 
these do not rely on uncontrolled approximations. How- 
ever, here are also some limitations. In order to study 
glassy behavior, we have to restrict the range of the con- 
trol parameters such that no crystallization effects occur. 
Furthermore, we need data with high numerical quality 
as input for MCT. For instance, the violation of the strict 
positive definiteness of the static correlation matrix St 
for a single k already leads to instabilities in the numer- 
ical solution of the MCT equations. Such problems typ- 
ically occur only for the atomistic computer simulation 
data at small k where both statistical noise and finite size 
effects may become dominant. Thus, we need a clearly 
defined way to smooth the numerical raw data such that 
S > is enforced. 

Now, let us come back to our specific model sys- 



3 



tem. We are able to reproduce available experimental 
data for Sk by performing MC simulations such that the 
deviations are smaller than the noise level of the ex- 
perimental data. Fig. Q] shows partial structure factors 
S°? for a binary mixture of dipolar particles in 2D at 
r = 110, 6 « 0.1, and x s = 0.48 measured experi- 



mentally by Ebert 01511 . Due to the finite polydisper- 
sity of both species of colloidal particles the value for 
6 is approximative, only. Also included are correspond- 
ing MC results which we have obtained as described 
in |Appendix A| with two modifications: First, we have 
produced these MC data by using only N = 600 par- 
ticles instead of N = 1600, and second, we have not 
smoothened the raw data. Thus, we observe that for low 
k (especially for S s k s ) the MC structure factors are su- 
perimposed by oscillations stemming from an interplay 
of uncertainties of h°P(f) at large r and the real space 
cutoff R for its Fourier transform. Besides these techni- 
cal artifacts, our MC results are in very good agreement 
with the experimental data. Note that there are no ad- 
justable parameters. For N = 600 particles, we were 
also able to reproduce experimentally measured partial 
structure factors within the parameter range 6 ~ 0.1, 
0.29 < x s < 0.48, and 51 < T < 281 which are not 
shown here. At smaller T, the influence of the hard-core 
interactions of the particles seems to become significant 
in the experimental data. For larger T, we would need 
larger system sizes and longer simulation runs in order 
to reproduce the experimental data. We can draw two 
important conclusions. First, forT > 51 the experimen- 
tal system of Konig et al. J2] is very well described by 
the idealized model potentials given by Eq. ([]]). Sec- 
ond, we can be confident that our MC simulations with 
1600 particles will yield sufficiently realistic results for 
our purposes, since the MCT glass transitions discussed 
below occur within the parameter range 80 < F < 160. 

Using integral equation theory, we have further an- 
alyzed several established choices for the bridge func- 
tions <ff(r) in Eq. The Percus-Yevick (PY) ansatz, 
which has successfully been applied to systems with 
hard core interactions, fails in describing MC simula- 
tion data for our system. For instance, for a binary mix- 
ture with 5 = 0.1, x s = 0.5, and F = 10, i.e. deep 
in the liquid regime, the PY closure underestimates the 
position of the main peak of S b k h by about 10 percent 
and overestimates the maximum value of this peak by 
almost 100 percent. For higher values of F we have 
not found stable numerical solutions for the PY struc- 
ture factors. Accordingly we do not know whether the 
quality of the PY theory for our model system becomes 
better or even worse. At least for F < 10, PY theory 
is much less reliable for dipolar potentials than for hard 



core interactions. The reason for this is not obvious. 
Another choice is the hypernetted chain (HNC) approx- 
imation d^(r) = 0. This simple ansatz describes the po- 
sitions of the main peaks of the partial structure factors 
correctly, but strongly underestimates their amplitudes 
(see below), as already noticed by Zahn 

G3l and Hoff- 
mann ilin . There are also more sophisticated closure re- 
lations. Such relations, however, typically include free 
adjustable parameters. This fact makes these methods 
not only computationally demanding, but also makes 
the optimization procedure unwieldy if more than one 
adjustable parameter is used. The Rogers-Young (RY) 
closure, for instance, which interpolates between PY 
and HNC, has in principle three fit parameters for the 
binary mixture. Furthermore, the quality of the corre- 
sponding results is not satisfactory. Fig. 3(a) in Ref. |Ql 
clearly shows that the RY closure underestimates the 
position of the main peak of the partial structure factors 
of the big dipoles. It also underestimates the amplitude 
of the oscillations beyond the main peak. 

For our binary dipole model in 2D, we have found 
empirically that using HNC with temperature T/2 in- 
stead of T leads to a systematic improvement in the de- 
scription of experimental and simulated data for S/< (see 
below). A similar observation was made earlier inde- 
pendently by Zahn ifioll . Even more, it was observed 
by Klapp |l2ll that also simulation results for specific 
dipolar systems in 3D can be well fitted by using HNC 
with T/2. These observations let us propose a simple 
empirical ansatz for the bridge functions: 



df(r) = - 



k R T ' 



(7) 



Since Eq. (Q) has no adjustable parameters, it is well 
suited for fast and systematic calculations. Eq. © will 
be called the T/2-HNC closure. 

In the following, we will use both, MC and T/2- 
HNC, for calculating the static structure input for MCT. 
Since the MC results for S 1 ? fit the corresponding ex- 
perimental ones very well, this ensures that the observed 
effects are predictions of MCT and not just artifacts of 
one of the specific techniques for the calculation of the 
static input. To calculate T/2-HNC structure factors nu- 
merically, we first eliminate the total correlations func- 
tions in Eqs. © and (0 by introducing the new func- 
tion y = h — C and apply then the Lado algorithm 113 ] 
with a real space cutoff R V" = 100 and N gr ;d = 4000 
grid points. The MC simulations are performed in the 
canonical ensemble with, if not explicitly stated other- 
wise, N = 1600 particles and 8 image boxes. Technical 
details and the procedure to smooth the raw data are ex- 



plained in Appendix A 
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Figure 2: (Color online) Partial structure factors for a binary mixture 
of dipolar particles in 2D at T = 120, S = 0.1, and x s = 0.5 from MC 
simulation with 1600 particles (solid lines) and for HNC (dashed lines 
in upper panel) and 772-HNC (dashed lines in lower panel). 



Fig. |2] shows MC simulated partial structure factors 
for the binary dipole model in 2D at T = 120, 6 = 0.1, 
and x„ = 0.5 compared to corresponding HNC and 772- 
HNC results. While the HNC closure underestimates 
the amplitude of the oscillations of the partial structure 
factors (see upper panel of Fig. [2]), the 772-HNC results 
are in a surprisingly good agreement with the MC data 
(see lower panel of Fig.[2j>. The agreement is not perfect 
since there are some smaller deviations, but the qualita- 
tive overall trends are even better than for the RY clo- 
sure. Remember that the RY closure underestimates the 
position of the main peak of the partial structure factors 
of the big dipoles and also underestimates the amplitude 
of the oscillations beyond the main peak (see Fig. 3(a) 
in Ref. H). The T/2-HNC data in Fig. [2] do not ex- 
hibit such qualitative drawbacks. Furthermore, we have 
found for all investigated examples within the control 
parameter range 0.1 < 6 < 0.6, 0.1 < x s < 0.9, and 
20 < F < 160 that, compared to the standard HNC clo- 
sure, the 772-HNC closure leads to a similar improve- 
ment in the description of MC data as shown in Fig. [2] 
In the next section, this trend will also be strongly sup- 
ported by the fact that using MC and 772-HNC structure 
factors as input for MCT will lead to compatible glass 
transition scenarios on very similar temperature scales. 

The MC data in Fig. [2] were obtained by simulating 
N = 1600 particles and smoothing the raw data as de- 
The positive definiteness S > 
However, there remain still some arti- 



scribed in Appendix A 



is guaranteed, 
facts at the lowest k values, see especially the kinks in 
S s k s for k < 0.8. By testing different extrapolation meth- 
ods we have found that these artifacts do not influence 
the glass transition scenarios shown below, and hence 



can be ignored. 

2.3. MCT equations 

The object of main interest for a statistical description 
of an isotropic and homogeneous classical fluid with 
M macroscopic components is the matrix O(f) of time- 
dependent partial autocorrelation functions of density 
fluctuations, <S>°f(f) (a, [3 = 1, . . . , M), at wave number 
k. The relation to static structure is given by 0(0) = 5, 
where the elements of the static structure factor matrix 
S obey lim^o, S'f = x a 5 a p. 

Starting point for the calculation of 0(f) is the exact 
Zwanzig-Mori equation. Assuming overdamped col- 
loidal dynamics, it reads 



T*(f)+S _1 ®(f) + 



Jo 



m(t-t')<b(t') = 0. (8) 



The positive definite matrix r of microscopic relaxation 
times will not be of further interest in this paper. MCT 
approximates m{t) by a symmetric bilinear functional 



m(t) = y-[*(r),0(r)]. 



(9) 



For a multicomponent liquid in d > 2 dimensions it 
reads explicitly 16J] 



(4. 
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[4k 2 p 2 - (k 2 + p 2 - ? 2)2](3-d)/4 ■ 



(12) 



c k denote the direct correlation functions already intro- 
duced above and = 27r^ 2 /T(d/2) is the surface of a 
unit sphere in d dimensions. F(x) is the gamma func- 
tion. 

The MCT Eqs. (l8ll-(fT2l can be solved numerically, 
only. For this we follow the standard procedure de- 
scribed in Ref. jf! and discretize the wave numbers. 
This leads to a finite, equally spaced grid of K points, 
k = (d d + k)Ak with % = 0, 1, . . ., K - 1 and < b d < 1. 
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The integrals in Eq. (fTTTt are then replaced by Riemann 



/-*oo pk+p 

I dp I dq. 

Jo J\k-p\ 



K-\ mm{K-\,k+p\ 

* (M)2 Z Z ■•■ (13) 

and Eq. ([8} represents a finite number of coupled equa- 
tions. For the offset we choose 02 = 0.303 J^l @|. 
The natural unit length is given by 1 / -\/n. The choice 
K = 250 and Ak = 0.2 turns out to be sufficiently accu- 
rate to avoid larger discretization effects. 

2.4. Glass transition lines 

To locate the glass transition point, we introduce 
the nonergodicity parameters (NEPs) given by F = 
lim^oo O(f). For our discretized model it can be proved 
1 14J] that F is (with respect to >) the maximum real sym- 
metric fixed point of 

I[X] =S-(S~ 1 +T[X,X]T 1 . (14) 

F can be calculated numerically by iterating Eq. (fT4l 
starting with X = S. 

For the binary dipole model with three independent 
control parameters, the glass transition takes place at 
a critical surface where F jumps from to F c > 0. 
This surface can be represented as r c {x s ,6). Fixing 
6 = 6 C defines a glass transition line (GTL). A general 
expression for its slope has been derived in Ref. |0]. 
Let (T c ,x c ^6 c ) be a critical point and (Ar, Ax s , A6) = 
(r - r c ,x s - x^,S - 6 C ). Then the following relation 
holds Hi 



dx s 



dcr/d(Ax s ) 



(15) 



(Ar,A.v,,,A<5)=0 



(x,„s)=(x-,6') do-/d(AT) 
where cr is the so-called separation parameter (!]]. 
Eq. <TT3T > can be evaluated in the so-called weak mix- 
ing limits x s — > and x/, — > by making use of the 
same perturbational approach as was performed for bi- 
nary mixtures of hard disks and hard spheres The 
few technical modifications which have to be taken into 
account for the binary dipole model are summarized 
Then Eq. (15[ serves as a powerful 
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tool for a fast and precise prediction of some qualitative 
properties of the GTLs. 

Let us also discuss some general properties of the 
GTLs. Our choice of T (Eq. (0) implies 

r c (x s , i) = r c (o,6) = r c (i,6) = r L (16) 

where FZ is the critical interaction parameter for the 
corresponding monodisperse system. Further, an inter- 
change of the roles of the big and small dipoles leads to 
the same physical scenario. Therefore it is 

r(x s ,5) = r(l-x s ,l/5). (17) 



Now, let us assume a small disparity in the susceptibil- 
ities, i.e. < (1 — <5) <K 1. Then, with Eq. ([Tol l we can 
write 

F(x s ,6) = T e Q + (l-6)T c 1 (xJ + Q.-Sfr 2 {x,) 
+0[(l-5f], (18) 

r7(0) = r^(i) = 0, ; = i,2. (19) 

Eq. ([PTT i implies for the coefficient functions 

Hfe) = -r^(i - Xs ), (20) 
r c 2 (x s ) = r 2 (i-x s )-r c 1 (i-x s ). (21) 

We can conclude: if rj(x s ) would be a non-vanishing 
function, then to leading order in (1 - S) we would ob- 
tain GTLs which, after subtracting Ty, are antisymmet- 
ric with respect to the equimolar composition x s = 1 /2. 
We could write (T c (x s , 6) - T C Q ) s - x s , 6) - F Q ) 

and r c (l/2, 6) = F c . If, however, there holds 



T\(x s ) = 



(22) 



for all x s while r c 2 (x s ) is a non-vanishing function, then 
Eqs. (TTHJ and d2T) imply 



r c (x s ,d) = r c (i-x s ,5) 



(23) 



which means that the GTLs become symmetric with 
respect to x s = 1/2 in the limit of small disparity in 
the susceptibilities. There holds also the (somewhat 
weaker) inverse statement. If Eq. ( l23l is assumed to be 
true for (1 - S) <k 1, then Eq. ( l20l implies Eq. ( l22l . For 
the slope s c (x s , 6) = (dT c /dx s )(x s , 6) of a GTL we can 
formulate further special (numerically testable) conclu- 
sions from Eq. d22"l) : 



(ds c /dd)(0, 1) = (ds c /dS)(l, 1) = 0. 



(24) 



Since the T/2-HNC closure relation has no free ad- 
justable parameters, their structure factors will be taken 
for the calculation of the slope d 15b . We will be able to 
produce numerical results with sufficiently high accu- 
racy to demonstrate that Eqs. ( 1231 and J241 I. and thus 
also Eq. ( 1221 ). are indeed valid for the binary dipole 
model. 

3. MCT results and discussion 

3.1. Glass transition lines 

Fig. [3] shows normalized slopes (dT c / dx a )(T c y l of 
the GTLs at x„ = 0, a = s,b, as functions of 1 1 6 calcu- 
lated by using Eq. (fT3T > and the T/2-HNC closure given 
by Eq. @. Because the case 5 — 1 is equivalent to a 
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Figure 3: (Color online) Normalized slopes of the GTLs for binary 
mixtures of dipolar particles in 2D at x s = and Xb = calculated 
by using Eq. i 1 5 1 and the 772-HNC closure given by Eq. {JJ. It is 
= 97.704 for the monodisperse limit. The inset presents the slopes 
on a much finer scale very close to 6 = 1 . 



one-component model, the slopes have to be zero at this 
point. Our numerical results clearly support this state- 
ment (see inset of Fig. |3). We observe that for < 8 < 1 
the slopes become always positive. This means that, 
compared to the corresponding one-component liquid, 
the presence of a small concentration of a second kind 
of particles always increases the critical interaction pa- 
rameter T c (x s , 8) at which vitrification sets in. For fixed 
temperature T and fixed total particle density n this im- 
plies that, compared to a given one-component system, 
the presence of a small concentration of a second kind 
of dipolar particles always increases the critical field 
B c and accordingly the critical value m c of the average 
magnetization per particle at which the system vitrifies. 
In this sense, mixing stabilizes the liquid state which is 
also called plasticization, especially in engineering lit- 
erature, see Ref. IU6I1 for example. 

Now, let us focus on the ^-dependence of the slope 
curve for x s = 0. Starting from zero at 8 — 1, the 
slope first increases upon decreasing 8. After exhibit- 
ing a maximum at 8+ = 0.444, the slope monotonically 
decays for asymptotically small 8. For Xb = the slope 
again starts from zero at 8 = 1 and increases strongly 
upon decreasing 8. As a further result, the inset in Fig. [3] 
demonstrates that the slope curves for x s = and Xb = 
become identical for (1 - 8) «; 1. This is already a first 
hint on the validity of Eq. d23b . i.e. on the symmetry 
of the GTLs with respect to i, = 1/2 for (1 - 8) «: 1. 
Furthermore, the numerical data indicate that the slopes, 
i.e. the derivatives of both slope curves vanish at 8 = 1, 
i.e. the data indicate that Eq. ( l24T i holds. 

The results in Fig.[3]already allow us to predict some 



Figure 4: (Color online) Glass transition lines for binary mixtures of 
dipolar particles in 2D with moderate susceptibility ratios calculated 
by using the 772-HNC closure, rjj is indicated by the dashed line. 



properties of the GTLs. Both, x s = and x s — 1, de- 
fine one-component models with the same critical in- 
teraction parameter r~. Hence, all GTLs must exhibit 
at least one maximum, since (<9F C /dx s )(x s = 0,(5) > 
and (dr/dx s )(x s = 1,8) = -(dT c /dx b )(x h = 0,8) < 
for < 8 < 1. Next, if for two values, 8\ and 82, it 
is (dr/dx s )(x s = 0,6i) < (dY c /dx s )(x s = 0,62) and 
(dr/dxbKxb = 0,6i) > (dr c /dx b )(x b = 0,62), i.e. 
(dT c /dx s )(x s = l,6i) < (dT c /dx s )(x s = 1,62), then 
the GTLs for 6\ and 82 have an odd number of inter- 
section points and therefore at least one crossing. Oth- 
erwise the number of intersection points is even. Since 
the slope curve s c (x s , 5) = (<9F C /dx s )(x s , S) at x s = de- 
creases monotonically with decreasing 5 for < 5 < 6+ 
while —s c {x s , 5) increases monotonically with decreas- 
ing 6 for all < 5 < 1 at x s = 1, two GTLs correspond- 
ing to 5 1 and 62 must have an odd number of crossings, 
if < 8\ < 82 < 8+. 

In the following, we present GTLs calculated both 
by using 772-HNC and MC structure factors as in- 
put for MCT Due to crystallization effects, the con- 
trol parameters for the MC data had to be restricted to 
0.1 < x s < 0.9 andO < 8 < 0.6. 

Fig.|4]shows GTLs for moderate susceptibility ratios 
8+ < 0.6 < 8 < 0.9 calculated by using T/2-HNC struc- 
ture factors as input for MCT. For < x s < 1, all GTLs 
are strictly above the monodisperse value F C Q and exhibit 
a single maximum. In addition, the GTLs do not in- 
tersect. The three lines for 0.7 < 8 < 0.9 are almost 
symmetric with respect to the equimolar composition 
x s = 1/2, i.e. the data strongly support Eq. d23l l and 
thus also Eq. (l22l . All these properties are compatible 
with the results for the slopes shown in Fig. [3] 

Fig- presents GTLs for 0.3 < 8 < 0.6 calculated 



7 




0.2 0.4 0.6 0.8 1 

x 



Figure 5: (Color online) Glass transition lines for binary mixtures of 
dipolar particles in 2D calculated by using 772-HNC structure factors 
(upper panel) and MC structure factors (lower panel). The dashed line 
in the upper panel indicates the 772-HNC value for r^. 



both by using 772-HNC and MC structure factors. The 
772-HNC curves again are strictly above the monodis- 
perse value T c for < x s < 1. For the MC curves, 
we can not test this statement strictly since we have 
no reasonable estimate for F£ due to crystallization ef- 
fects. Therefore, we focus first on the 772-HNC results. 
While the GTL for 5 = 0.6 exhibits only a single rather 
fiat maximum, the other GTLs with 0.3 < 5 < 0.5 
develop two local maxima with a local minimum at 
x s ~ 0.6 in between. Due to the occurrence of this lo- 
cal minimum for 5 = 0.5 and the increase of both slope 
curves s c (x s ,6) at x s — and -s c (x s , 5) at x s — 1 with 
decreasing 5 for 5+ < 5 < 1, the curves for 5 = 0.5 and 
5 = 0.6 have two intersection points within the interval 
< x s < 1 . If, however, we plot a pair of the shown 
GTLs with 0.3 < 5\ < 62 < 0.5, then we obtain exactly 
one intersection point within the interval < x s < I. 
For the pair 5i = 0.3 and 62 = 0.4 the existence of 
an odd number of intersection points is enforced by the 
facts \(dr c /dx s )(x s = 0,501 < \(dT c /dx s )(x s = 0,52)| 
md\(dr/dx s )(x s = l,5i)| < |(3T7d*,X*i = 1,^)1 
following both from the GTLs by numerical differenti- 
ation and also from the slopes in Fig. [3] since 5, < 5+ 
for i = 1,2. The upper panel of Fig. [5] shows that only 
one intersection point occurs for the two GTLs. Now, 
let us consider the MC results in Fig. [5] The data are 
numerically less precise and the occurring extrema in 
the curves are not as clearly pronounced as for the 772- 
HNC results. Furthermore, the MC curves do not ex- 
hibit the same number of intersection points as the cor- 
responding 772-HNC ones. Nevertheless, the MC data 
clearly support at least the qualitative ^-dependencies 
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Figure 6: (Color online) Glass transition lines for binary mixtures of 
dipolar particles in 2D with larger disparities in the susceptibilities 
calculated by using 772-HNC structure factors (upper panel) and MC 
structure factors (lower panel). The dashed line in the upper panel 
indicates the 772-HNC value for T C Q . In the neighborhood of x s = 0.8 
and S = 0. 1 we had problems in finding stable numerical solutions for 
the r/2-HNC structure factors, and thus there are no data points for 
r c in this region in the upper panel. 

of the GTLs predicted by the 772-HNC closure. Note 
that the scales for T c and its total variation are also sim- 
ilar in both approaches. 

Fig. [6] shows GTLs for larger disparities in the sus- 
ceptibilities, 0.1 < 5 < 0.3, calculated both by using 
772-HNC and MC structure factors as input for MCT 
Let us first focus on the 772-HNC data. Again, all 
GTLs are strictly above the monodisperse value r c Q for 
< x s < 1. The GTL for 5 = 0.3, which is also 
shown in Fig.|5]on a finer scale for T c , exhibits two lo- 
cal maxima at x s * 0.25 and x s w 0.88 and in between 
a clearly pronounced local minimum at x s ~ 0.6 (we 
have obtained these numbers via spline-interpolations). 
Upon decreasing 5, the amplitude of the maximum at 
x s « 0.88 strongly increases while the maximum at 
x s « 0.25 and thus also the minimum at x s ~ 0.6 be- 
come less pronounced. As a result, the GTL for 5 = 0.2 
exhibits only a single maximum at 5 ~ 0.89. The GTL 
for 5 = 0.1 seems to be qualitatively similar, but in the 
neighborhood of x s = 0.8 and 5 = 0.1 we had problems 
in finding stable numerical solutions for the 772-HNC 
structure factors, and thus we have no data points for F c 
in this region. As already predicted by the slopes pre- 
sented in Fig. [3] each pair of the shown GTLs exhibits 
an intersection point within the interval < x s < 1 . The 
MC results in Fig.|6]clearly support the trends predicted 
by the 772-HNC closure. Note especially the similarity 
of the intersection points of the curves. Further, both 
approaches yield almost identical scales for F c and its 
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Figure 7: (Color online) Self-intermediate scattering functions of 
Konig et al. Ql measured experimentally at 6 ss 0.1 and x s as 0.3 
(upper panel) and x s a 0.5 (lower panel). The functions were av- 
eraged over big and small particles. k max denotes the wave number 
corresponding to the first maximum of g bb (r). SeeRef. Q] for details. 
The dotted horizontal lines mark the value 0.6 (see text). 



total variation. 

The common feature of all T/2-HNC results in 
Figs- EJ-EJis the fact that there holds strictly 



T c {x s ,5)>Ti 



for all tested < 6 < 1 and < x s < 1 . The correspond- 
ing MC results in Figs.|5]and|6]are also consistent with 
these findings. Thus we can conclude: for the binary 
dipole model in 2D, MCT predicts always a stabilization 
of the liquid state due to mixing. Further, the amplitude 
of this plasticization effect is quite large. For instance, 
the variation of F c in Fig.[5]is of the order of 10 percent. 
The GTLs for 6 = 0.1 (Fig. [6) vary by even more than 
50 percent. Tests of MCT in the last two decades have 
shown that for model systems with smooth pair poten- 
tials the critical temperatures predicted by MCT differ 
from those determined from computer simulations by 
about a factor of two, i.e. 100 percent |l|]. This does not 
necessarily imply for the binary dipole model that the 50 
percent variation of F c for 6 = 0.1 can not be observed 
experimentally. If T c exp (x s ,6 = 0.1) ~ 2Y c (x s ,6 = 0.1) 
for < x s < 1, then the experimentally determined crit- 
ical vale T c exp would vary by about 50 percent, as well. 
Thus, the predicted plasticization effect should be ex- 
perimentally testable. 

3.2. Comparison to experiments 

Of course, it would be interesting to determine 
T c (x s , 6) experimentally, which has not been done so far. 
Another possibility to check our prediction, as already 
discussed for binary mixtures of hard disks and hard 



spheres in Refs. flO], is the strong variation in the a- 
relaxation time scale (T rel )°f upon changing x s or 6. Let 
us fix 6 — 0.1 (which is approximately the experimen- 
tal value J2l) and choose two compositions x { P — 0.3 



and 



0.5. Now, let us further fix F such that 



< [Fixf^) 



totic scaling law (T rel ') c f 



H/r «. 1 for i = 1, 2. With the asyrnp- 

([r c (x s ,6)-r]/ry^"V 



and the fact T c (x ( s n ,6) < T c (x ( P,6) (see Fig.© we ob- 



tain (T re/ )f(r,x< 2) ,<5) 



« (^ l )f(T,x s l) ,6) 



Thus, we 
0.1 that binary 



expect for this chosen value of F and 6 
mixtures with x s - xf^ = 0.5 will behave more liquid- 
like than corresponding mixtures with x s — x^ — 0.3, 
i.e. their structural relaxation is faster. 

Next we consider the experimental results of Konig 
et al. J2l for the self-intermediate density correlators 
shown in Fig. Q which were obtained for 6 as 0.1 and 
two different compositions x s ~ 0.3 and x s as 0.5. Note 
that the experimental values for x s of Konig et al. J2l are 



not as accurate as the corresponding ones of Ebert 115 ] 
(see caption of our Fig.[TJ. From Fig. [6] we find F c (x s = 
03,6 = 0.1) = 104 < r c (x s = 0.5,(5 = 0.1) = 115 in 
case of the static input from our MC simulation. On the 
other hand, for the structural relaxation time r(r, x s , 6) 

self 

defined by <EC' J (t, T, x s , 6) = 0.6 we can read off from 



(25) Fig. □ that T exp (T 



exp\ 

Texp(X = 200, X s * 

t(Ti,x s ,S) < t(T 2 ,x 
be 



226, x s ~ 0.5, S ~ 0.1) = 569.? « 
0.3,6 ~ 0.1) = 14513*. Since 
t, 6) for Fi < F 2 < T c (x s , 6), it must 



T exp (T = 200, x s * 0.5, 6 * 0.1) « T exp (T = 

200, x s * 0.3,6 -0.1). (26) 

Let r c etp (x s ,6) be the critical value at which an ex- 
trapolation of T exp (T,x s , 6) diverges as a power law 
T exp (T,x s ,6) ~ [r c exp (x s ,6)-T]-^ s \ Assuming that 
the exponent and the proportionality factor of that power 
law do not change sensitively under a change of x s as 
0.3 to x s ~ 0.5 (for fixed 6 ~ 0.1), then the in- 
equality (l26b implies that the distance of T = 200 to 
T° exp (x s as 0.3,6 as 0.1) is smaller than the distance 
as 0.5,6 as 0.1). Consequently, we have 
0.3,6 as 0.1) < r c exp (x s as 0.5,6 as 0.1), in 
agreement with the MCT predictions above. This in- 
equality between the experimental values for T c exp is also 
consistent with the fact that the mixture with x„ as 0.3 
and T = 200 already shows a plateau whereas the mix- 
ture with x s as 0.5 and F = 226 does not. 

This first attempt to compare our results with the ex- 
perimental ones indicates that MCT is able to predict at 
least qualitatively the change in the relaxation behavior 
upon composition change from x s as 0.3 to x s as 0.5. Of 
course, it is necessary to strongly extend the parameter 



to T c exp (x, 
F fV p(x s a 
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range for x s , but also for 5, to test the validity of our 
MCT predictions more systematically. Performing such 
a systematic test has also to account for the overestima- 
tion of the temperature scale at which vitrification sets 
in. For 8 = 0.1 and x s = 0.3 MCT yields P * 105 by 
both using 772-HNC and MC structure factors as input. 
In the upper panel of Fig. [7] however, the development 
of a plateau as a precursor of a possible glass transition 
only becomes clearly visible in the experimental data 
for r > 200. Hence, MCT overestimates the scale for 
the glass transition temperature by about a factor of two, 
as well-known for other systems. 

4. Summary and conclusions 

The major motivation of our study has been to ex- 
plore the effect of mixing on the ideal glass transition 
of a binary mixture of dipolar colloids in 2D and the 
comparison with corresponding experimental data. The 
calculations were done in the framework of MCT. 

The MCT equations were solved by using the static 
structure factors from MC simulations and integral 
equation theory. For the latter we have used the T/2- 
HNC closure relation because it does not involve ad- 
justable parameters and fits the MC results for the par- 
tial static structure factors rather well. The very good 
agreement of the partial static structure factors from ex- 
periment J2l with those from our MC simulation using 
dipolar point particles has demonstrated that the dipo- 
lar colloidal system is well modeled by point particles 
in 2D interacting only via the induced dipole potentials 
([TJ, at least in the parameter range where the glass tran- 
sition occurs. 

The main result we have found is that the effect of 
mixing on the glass transition depends strongly on the 
liquid system. In contrast to binary mixtures of hard 
disks and hard spheres |@,01, MCT for the binary liq- 
uid of point dipoles predicts always a stabilization of 
the liquid state due to mixing. More explicitly, for fixed 
temperature T and fixed total particle density n the crit- 
ical magnetic field B c and accordingly the critical value 
m c of the magnetization per particle at which the binary 
dipolar liquid vitrifies is always larger, compared to a 
one-component system. The amplitude of this predicted 
plasticization effect is quite large, the GTLs for 8 — 0.1 
shown in Fig. [6] vary by more than 50 percent. 

The ff-relaxation times of the experimentally ob- 
tained self-intermediate density correlators shown in 
Fig.|7]for x s ~ 0.3 and x s ~ 0.5 clearly support the pre- 
dicted plasticization effect. This hints that MCT may 
be able to predict qualitatively the change in the relax- 
ation behavior upon composition changes. The exper- 
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Figure 8: (Color online) Partial structure factors Sf a for binary mix- 
tures of dipolar particles in 2D at T = 95 and x s = 0.5 calculated by 
using the 772-HNC closure. 

imental data also show that MCT systematically over- 
estimates the scale for the glass transition temperature 
by more than a factor of two, which is a well-known 
phenomenon of MCT (l|] . 

At this point, let us also speculate about the physical 
origin of the predicted plasticization effect. As pointed 
out by Hoffmann et al. 0], the smaller particles in 
the binary dipole model tend to form clusters whereas 
the bigger particles remain more homogeneously dis- 
tributed in space. This partial clustering manifests itself 
by the occurrence of a prepeak in the partial structure 
factors S k of the smaller particles at low k. Further- 
more, the cited clustering effect seems to occur as a 
generic feature of the binary dipole model: all partial 
structure factors SI' with < 8 < 1 and < x s < 1 
used for calculating the glass transition lines in Figs.@} 
exhibit a prepeak. On the other hand, MCT predicts 
plasticization for all < 8 < 1 and < x s < 1 where the 
only model specific input to MCT are the partial struc- 
ture factors 5f . Thus, it could be that the partial clus- 
tering of the smaller particles is the physical origin for 
the plasticization effect predicted by MCT. 

To investigate this conjecture in more detail, we have 
calculated 772-HNC structure factors S'f for F = 95, 
x s = 0.5, and 8 = 1.0, 0.8, 0.6, and 0.4 from which the 
bb and ss elements are shown in Fig. [8] The case £=1.0 
is a monodisperse system close to its MCT glass transi- 
tion point which shows no clustering and by definition 
satisfies S b k — S s k s . Now, by decreasing 8 we observe 
five qualitative effects on S? a : (i) a growing prepeak in 
S I s at 2 < k < 5, (ii) a decrease in for < k < 3, (iii) 
a decrease in the amplitude in the main- and the subse- 
quent peaks in both S b k h and Sl s , (iv) a shrinking of the 
k scale for the oscillations of S h k , and (v) a stretching of 
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the k scale for the oscillations of Sf s . The plasticization 
effect predicted by MCT must be a result from an inter- 
play of these five effects. Since these effects are com- 
peting, no definite conclusion can be drawn concerning 
the plasticization effect. Thus, at least from our analy- 
sis here it is not obvious whether there is a connection 
between clustering and plasticization or not. 

We have also discussed some general properties of 
the glass transition diagram and its GTLs. For small 
disparities in the susceptibilities, i.e. for 5 close to one, 
the GTLs are symmetric like for binary hard disks SQl . 
Two GTLs with d\ and 82 below 6+ ~ 0.444 must have 
an odd number of intersection points. This prediction is 
only based on the application of the analytical expres- 
sion for the slope (dT° /dx s )(x s , 8) at x s = and x s = 1. 
Furthermore, depending on 5 the GTLs exhibit either 
one maximum, for 6 closer to one or 6 very small, and 
two maxima with a minimum in between, for intermedi- 
ate values of 5. This behavior also differs qualitatively 
from that of binary hard disks in 2D and binary hard 
spheres in 3D. 

As a side result to be mentioned, with the T/2-HNC 
closure we have suggested a new bridge function ap- 
proximation which is useful for calculating the static 
structure functions of dipolar systems. Of course, our 
finding is only of empirical nature. 

The various properties of the GTLs of binary liquids 
of point dipoles in 2D predicted here in the framework 
of MCT are a challenge to be tested in more detailed 
numerical simulations and by extending the present ex- 
perimental approach. But the theoretical approach can 
be extended as well. For instance, for the present model 
it will be interesting to study the fully time-dependent 
density correlators. A further interesting topic will be 
the discussion of the glass transition behavior of bi- 
nary hard disks with dipolar interactions at intermediate 
and high packing fractions, where the mixing effects for 
dipolar point particles and those for hard disks |@] may 
interfere. Such investigations are work in progress and 
will be published elsewhere. 
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Appendix A. Monte Carlo simulation 

We simulate N = 1600 particles within a cubic box 
of dimensionless length L = y/N - 40 with the stan- 
dard periodic boundary conditions for particle motion. 
This corresponds to setting the total particle number 
density to 1 since the magnetic field B and the suscep- 
tibility ratio 6 are the only physical control parameters, 
see Eq. (0. We arrange 8 virtual copies of the simula- 
tion box around it. These so-called image boxes con- 
tain at every simulation step exactly the same particle 
configurations as the main simulation box. To calcu- 
late the potential energy of a particle inside the main 
simulation box, the interactions with all particles within 
the main and the 8 image boxes are taken into account. 
This is necessary due to the long -ranged 1/r 5 pair po- 
tentials. However, these potentials are absolutely in- 
tegrable in 2D and thus we can calculate the potential 
energies without advanced techniques like Ewald sum- 
mation. A trial move for a chosen particle consists of a 
local random displacement Ax, Ay € [-£, s], s = 0.15 
of the considered particle (and also of its copies in the 
image boxes) followed by calculating the potential en- 
ergy cost Am for this displacement and choosing a ran- 
dom number z e [0,1]. The trial move is accepted if 
z < e\p(-Au/ksT), and rejected, otherwise. Typical 
acceptance probabilities are between 0.3 and 0.5. A so- 
called sweep is a sequence consisting of exactly one trial 
move for every particle in the main simulation box. 

The calculation of the static correlation matrix S* 
consists of four steps. First, we randomly distribute 
the N = 1600 particles. Second, we perform 10000 
sweeps for equilibrating the system. As third step, we 
perform 20000 sweeps to record a histogram for the ra- 
dial distribution functions g°Hr). For this purpose, r is 
discretized as a Lado grid 111311 with a real space cutoff 
R = L/2 = 20 and N gr id = 500 grid points. By calculat- 
ing the Fourier transform of h a ^(r) = g a ^(r) - 1 accord- 
ing to the Lado method H 1 3fl we obtain the raw data for 
h. As last step, we smooth the data. The raw data for 
h°f plotted as functions of k are superimposed by os- 
cillations stemming from an interplay of uncertainties 
of h a P(r) at large r and the real space cutoff R. These 
oscillations are most significant at small k. We smooth 
them out by taking the arithmetic average of the upper 
and lower envelope functions. With this result we cal- 
culate the raw data for Now, for the first few of the 
lowest k grid points the positive definiteness of may 
be violated. We repair this as follows: let / be the small- 
est index for which S*. > holds for all j > I. We set 
Sk, = Sk M for / < /. The whole procedure can be fully 
automatized. 
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Appendix B. Weak mixing limit 

The evaluation of Eq. (Q3J) at x s = and jq, = 
follows exactly the steps described in the Appendix of 
Ref. |@] with the following modifications: the pack- 
ing fraction (p has to be replaced by the variable T. 
Eqs. (A29) and (A31) in Ref. 0] have to be replaced 
by the zeroth order T/2-HNC closure 

ln[(/i (( W) +1] = (fc (0) ) ff/ V) - (c (0) ) ff/ V) 

-(uf} f ) ap (r), (B.l) 



r „(0) .a/3, , _ XaXfi T 

<- U eff> (0 - { ^(0)}3/2 x 2 73' 



and the first order T/2-HNC closure 

(A (1) )^(r) = {(h m ) afl (r)+ l}{(h m )^(r) 

-(c (1) rV) - (u%) aP (r)}, (B.3) 

(H^)^(r) = 2{1 - S}(uf} f )^(r). (B.4) 

Furthermore, Eq. (A32) in Ref. JfJ has to be replaced 
by n (0) = 1 and n (1) = 0, since we have chosen 1/ -\/n as 
unit length. 
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